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ABSTRACT 


The  contribution  of  inertia  terms  towards  load  bearing  capacity, 

the  effect  of  convection  terms  on  the  temperature  distribution  and  the 

temperature  variation  in  high  and  low  viscosity  lubricants,  are  considered 

for  a  plane  slider  bearing.  The  fluid  was  to  be  assumed  incompressible 

of  constant  but  arbitary  viscosity.  Through  a  stream  function  approach 

and  a  power  series  expansion,  the  momentum  equation  is  converted  into  a 

set  of  ordinary  differential  equations  which  are  numerically  solved.  The 

pressure  distribution,  load  capacity,  shear  stress  and  the  increase  in 

* 

load  capacity  is  calculated  for  various  modified  Reynolds  numbers,  Re. 

The  energy  equation  is  also  solved  numerically  with  and  without  convective 

terms.  The  temperature  distribution  is  obtained.  It  was  found  that  the 

presence  of  the  inertia  terms  increases  the  load  capacity.  This  effect 

* 

becomes  important  at  high  Re.  The  analysis  indicates  that  the  convection 
terms  can  not  be  neglected  and  that  the  lubricating  oils  of  high  vis¬ 
cosity  develop  larger  temperature  variation  in  the  film  than  those 
of  low  viscosity. 
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CHAPTER  I 

STATEMENT  OF  THE  PROBLEM  AND  REVIEW  OF  RELEVANT  LITERATURE 

1.1  INTRODUCTION 

The  classical  theory  of  hydrodynamic  lubrication  is  based  upon 
the  assumption  that  the  dominating  influence  is  that  of  lubricant  viscosity 
In  the  narrow  space  between  two  inclined  surfaces  in  relative  motion,  the 
pressure  generated  in  a  fluid  film  by  viscous  forces  tend  to  keep  the 
surfaces  apart,  thus  forming  a  bearing.  The  most  favourable  combinations 
of  geometry  and  kinematics  of  such  system  have  been  much  sought  after. 
Predictions  are  usually  required  of  the  main  dependent  variables  pressure, 
load  capacity,  shear  stress,  mass  flow  rate  and  temperature  distribution; 
these  are  influenced  by  the  geometry  of  the  system,  the  kinematics  of 
the  system  and  the  properties  of  the  lubricant. 

Specifically,  the  basic  differential  equationsfor  the  study  of 
hydrodynamic  lubrication  are  the  continuity  equation,  the  momentum 
equations,  the  energy  equations  and  the  equation  of  state.  These  equations 
interdependent  as  they  are  through  the  physical  properties  of  the  lubri¬ 
cant,  are  extremely  complicated  mathematically.  Several  simplifying  as¬ 
sumptions  have  to  be  made  before  a  solution  could  be  obtained. 

The  hydrodynamic  lubrication  may  be  broadly  classified  under 
two  groups,  depending  on  the  type  of  the  lubricant,  namely 
(a)  compressible  hydrodynamic  lubrication 
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(b)  imcompressible  hydrodynamic  lubrication. 

Each  of  the  above  two  groups  could  be  further  subdivided  into 
(i )  self  acting  bearing 
(ii)  externally  pressurized  bearing. 

In  a  bearing  different  surface  configurations  are  possible  and  prevalent, 
but  the  two  most  common  configurations  are 
(i)  the  journal  bearing 
(ii)  the  thrust  bearing. 

In  this  thesis,  the  analysis  will  be  confined  to  incompressible 
hydrodynamic  lubrication  and  the  assumed  configuration  will  be  a  plane 
slider,  representing  a  thrust  bearing. 

1.2  INERTIA  EFFECTS  AND  TEMPERATURE  VARIATION  IN  SLIDER  BEARING 

The  Reynolds  equation  [1]  is  based  on  the  assumption  that  the 
inertia  forces  are  negligible  compared  to  the  viscous  forces.  The  im¬ 
portance  of  the  inertia  terms  relative  to  the  viscous  terms  in  the  Navier 

Stokes  equations  can  be  characterized  by  a  dimensionless  parameter  re- 

* 

f erred  to  as  the  modified  Reynolds  number,  Re. 

* 

In  the  operation  of  most  of  the  bearings  the  Re  is  very  small, 

therefore,  the  inertia  effects  could  be  easily  ignored.  However,  for 

bearings  operating  at  high  speeds,  particularly  if  the  lubricant  is  of 

*• 

low  kinematic  viscosity,  the  Re  may  assume  values  near  or  exceeding  unity 
Therefore,  the  inertia  forces  become  comparable  with  the  viscous  forces 
and  the  inertia  effects  may  no  longer  be  negligible. 

BRAND  [2]  has  shown  that  if  the  Reynolds  number  pRcohm/y  is  of 
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the  same  order  as  R/h  ,  then  the  inertia  forces  will  be  of  the  same  order 
as  the  viscous  forces.  However,  if  the  Reynolds  number  becomes  too  large, 
turbulence  will  develop  in  the  flow  and  the  governing  equations  and  the 
analysis  will  have  to  be  modified  suitably. 

Another  assumption  made  in  the  classical  hydrodynamic  lubri¬ 
cation  is  that  the  effects  of  temperature  variation  in  the  lubricant  film 
are  negligible.  Therefore,  the  solution  to  the  momentum  equation  is  ob¬ 
tained  for  constant  viscosity.  However,  in  reality,  energy  is  generated 
by  viscous  dissipation  which  results  in  the  heating  of  the  lubricant  in 
the  bearing,  therefore  decreasing  the  magnitude  of  the  lubricant  viscosity 
Hence  the  solution  obtained  from  the  momentum  equation  for  constant  vis¬ 
cosity  will  be  in  error.  In  order  to  obtain  an  approximate  solution  in 
accordance  with  the  actual  operation  of  the  bearing  some  mean  viscosity 
value  was  employed  in  the  solution  of  the  momentum  equation.  The  pre¬ 
diction  of  the  mean  viscosity  value  was  made  from  the  inlet  and  outlet 
temperature  conditions. 

Several  authors  have  attempted  to  include  the  viscosity  vari¬ 
ation  into  account  by  allowing  the  viscosity  variation  in  x-di recti  on 
only.  This  results  in  a  much  simpler  form  of  the  momentum  and  the  energy 
equations.  These  differential  equations  were  solved  simultaneously  for 
the  imposed  boundary  conditions. 

The  importance  of  the  lateral  temperature  variation  depends 
strongly  on  the  thermal  boundary  conditions  imposed  on  the  film  boundaries 
If  both  boundaries  are  taken  as  adiabatic,  all  the  energy  generated  by 
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viscous  dissipation  would  have  to  be  carried  away  by  the  fluid.  If  both 
boundaries  were  maintained  at  constant  temperature,  greater  lateral 
temperature  gradient  would  develop  and  heat  transfer  across  both  boundaries 
would  occur. 

1.3  REVIEW  OF  RELEVANT  LITERATURE 

The  retention  of  inertia  terms  in  the  simplified  Navier-Stokes 
equations  for  the  slider  bearing  yields  a  non-linear  differential  equation. 
The  influence  of  these  non-linear  inertia  terms  have  been  investigated  by 
a  number  of  authors  adopting  approximate  methods. 

SLEZKIN  and  TRAG  [3]  have  proposed  that  the  inertia  terms  are 
averaged  across  the  film  thickness.  Therefore  the  inertia  terms  become 
a  function  of  longitudinal  direction  only,  and  the  resulting  equation 
may  be  readily  integrated. 

Several  authors  [4,5]  have  employed  this  technique  to  obtain 
the  solution  for  various  bearing  configurations.  The  averaged  inertia 
approximation  results  in  a  parabolic  velocity  profile  which  need  not  al¬ 
ways  be  applicable. 

KAHLERT  [6]  developed  a  method  of  solving  the  differential 
equation  in  which  it  is  assumed  initially  that  the  inertia  forces  are 
small  in  comparison  with  the  viscous  forces.  Hence,  their  effects  may 
be  calculated  as  a  small  correction  to  the  expression  obtained  from 
purely  viscous  consideration.  The  noticed  difference  between  the  re¬ 
sults  of  Kahlert  and  the  other  authors  can  be  traced  to  the  approximation 
involved  in  assuming  the  inertia  forces  as  a  small  fraction  of  the  viscous 


' 
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forces . 

MILNE  [7]  has  considered  a  solution  to  the  complete  equation 
of  motion  in  cylindrical  co-ordinates  without  using  the  thin  film  approxi¬ 
mation.  His  technique  consists  of  a  series  expansion  of  the  stream 
function  with  reciprocal  of  kinematic  viscosity  as  the  expansion  para¬ 
meter.  The  zero  order  term  in  the  series  corresponds  to  the  inertia¬ 
less  case  since  1/v  -+  0  for  negligible  inertia.  In  a  later  review  he 
points  out  that  this  procedure  is  identical  to  that  developed  by  Kahlert. 

SNYDER  [8]  has  obtained  a  more  exact  solution  by  taking  into 
consideration  the  variation  of  the  inertia  effects  across  the  film  as 
well  as  in  the  direction  of  the  flow.  In  this  method  the  stream  function 
approach  is  introduced.  The  stream  function  is  expressed  as  a  power 
series  in  6  which  is  the  function  of  h(x),  whereas  the  coefficients  are 
assumed  to  be  a  function  of  n  which  in  turn  depends  on  y  co-ordinate. 

The  coefficients  are  obtained  by  solving  a  set  of  differential  equations. 

In  presenting  the  solution  SNYDER  solved  for  only  the  first  two  terms  of 
the  power  series.  The  treatment  presented  here  using  the  numerical 
method  indicates  that  it  is  necessary  to  include  at  least  the  first  four 
terms  to  obtain  more  accurate  results. 

WOODHEAD  and  KETTLEBOROUGH  [9]  have  obtained  the  solution  of 
the  differential  equation  by  matrix  methods.  They  have  introduced  the 
concept  of  field  matrix  which  reduces  the  problem  of  determining  the 
three  field  matrices  [U],  [V]  and  [P],  The  solution  is  obtained  by 
iterative  procedure  until  the  boundary  conditions  are  satisfied.  Generally, 


•  29H92  *T9W0<J  9riJ 


the  solution  by  matrix  methods  is  long  and  cumbersome  and  takes  sub¬ 
stantial  computer  time. 
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HUNTER  and  ZIENKIEWICZ  [10]  considered  the  momentum  and  the 
energy  equation  in  the  coupled  form  for  the  isothermal,  as  well  as  the 
adiabatic  bearing  surfaces.  The  viscosity  of  the  lubricant  is  regarded 
as  a  function  of  temperature.  In  their  analysis  they  have  omitted  the 
inertia  terms  in  the  momentum  equation  and  the  v  8T/3y  term  in  the  energy 
equation.  The  solution  to  the  differential  equations  is  obtained  by  a 
finite  difference  method  using  an  iterative  procedure.  The  omission 
of  inertia  terms  indicates  that  the  solution  is  valid  for  very  low 
modified  Reynolds  numbers. 

The  latest  investigation  was  carried  out  by  HAHN  and  KETTLE- 
BOROUGH  [11].  The  coupled  form  of  the  momentum  and  the  energy  equation 
is  considered  for  the  isothermal  surfaces.  The  viscosity  of  the  lubricant 
is  regarded  as  a  function  of  temperature  and  pressure.  They  have  re¬ 
tained  the  inertia  terms  in  the  momentum  equation  and  the  convective 
terms  in  the  energy  equation.  The  solution  is  obtained  through  matrix  methods 
similar  to  the  procedure  mentioned  in  [9]. 

On  examination  of  the  literature  pertinent  to  the  incompressible 
hydrodynamic  lubrication  it  seems  apparent  that  the  complete  solution  to 
the  momentum  and  the  energy  equations  is  difficult  due  to  their  non¬ 
linearity  and  coupled  form.  To  achieve  this  object  it  seems  necessary  to 
introduce  some  form  of  transformation  that  would  simplify  the  differential 
equations  so  that  a  simpler  numerical  solution  is  possible.  It  appears 


that  stream  function  methods  have  great  potential  in  the  study  of  lubri¬ 
cation  films.  It  is  felt  that  approach  on  these  lines  would  further  en¬ 
hance  the  knowledge  of  film  lubrication. 

1.4  STATEMENT  OF  THE  PROBLEM 

The  purpose  of  the  present  investigation  is: 

(a)  to  investigate  the  contribution  of  inertia  forces  for  different 
modified  Reynolds  numbers, 

(b)  to  consider  the  effect  of  convective  terms  in  the  energy 
equation,  and 

(c)  to  consider  the  temperature  variation  in  the  lubricating 
film  of  high  and  low  viscosity  oils. 

The  bearing  configuration  considered  is  a  plane  slider. 

In  obtaining  the  solution,  the  assumption  of  constant  viscosity 
will  be  made.  This  results  in  uncoupling  of  the  momentum  equation  from 
the  energy  equation  so  that  an  independent  solution  for  the  velocity 
distribution  can  be  obtained. 

The  uncoupled  momentum  equation  is  then  transformed  through  a 
stream  function,  yielding  a  set  of  ordinary  differential  equations.  These 
equations  are  then  solved  numerically.  The  dependent  variables  such  as 
pressure  distribution,  load  capacity,  and  the  shear  stress  at  the  moving 
surface  are  computed  for  various  values  of  the  modified  Reynolds  number, 

'jlf 

(0.1  £  Re  £  1.0).  The  increase  in  load  capacity  due  to  inertia  effects 

★ 

is  then  computed  for  the  above  Re. 

The  solution  to  the  energy  equation  is  also  obtained  numerically 


8 


with  and  without  convective  terms.  The  temperature  variation  in  the 
film  and  the  temperature  gradient  at  the  moving  and  the  stationary 
surfaces  are  calculated.  In  addition,  these  results  have  also  been  ob¬ 
tained  for  bearings  using  high  and  low  viscosity  oils.  For  these  cases 
the  energy  equation  considered  included  the  convective  terms. 

The  relevant  results  have  been  discussed  in  the  final  chapter. 


. 


9 


CHAPTER  II 

THE  GOVERNING  EQUATIONS 

2.1  GENERAL  FORM  OF  THE  EQUATIONS 

Consider  the  Navier-Stokes  equations  and  the  energy  equation 
[12]  for  the  slider  bearing  geometry  of  Fig.  1.1. 


FIG.  1.1  SLIDER  BEARING  GEOMETRY 


The  Navier-Stokes  equations  and  the  energy  equation  are  simpli¬ 
fied  under  the  following  assumptions: 

1.  The  problem  is  two  dimensional. 

2.  Steady  state  has  been  reached. 

3.  The  effects  of  thermal  and  elastic  distortions  are  neglected. 
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4.  The  lubricant  density  and  viscosity  remains  constant. 

5.  The  thermal  conductivity  and  specific  heat  of  the  lubricant 
are  constant. 

6.  The  flow  is  laminar. 

7.  Cavitation  effects  may  be  ignored, 

8.  Lubricating  film  is  very  thin,  h/L  «  1, 

9.  Order  of  magnitude  considerations  may  be  used  to  simplify 
the  Navier-Stokes  and  the  energy  equations. 

10.  The  inlet  temperature  of  the  lubricant  may  be  taken 


constant. 


The  governing  equations  may  be  written  as: 


Momentum  Equations 


(2.1) 


(2.2) 


Continuity  Equation 


(2.3) 


Energy  Equation 


(2.4) 
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The  boundary  conditions 


u  =  U,  v  =  0  at  y  =  0 

u  =  0,  v  =  0  at  y  =  h(x) 

\ 

p  =  0  at  x  =  0 

> 

p  =  0  at  x  =  L 

/ 

T(0,y)  =  Ts 
T(x,0)  =  Ts  \ 

T(x,h)  =  T$  ) 


(2.5) 


(2.6) 


(2.7) 


2.2  EQUATIONS  IN  DIMENSIONLESS  FORM 

The  governing  equations  may  also  be  written  in  the  dimension¬ 
less  form  by  introducing  the  following: 

x  =  xL 

y  =  yh 

u  =  u  U 

v  =  vV 


' 


■ 


V 


\ 
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p  =  p  pU 


2 


The  equations  (2.1)  through  (2.7)  can  be  written  in  the  fol¬ 
lowing  dimensionless  form: 


.  JE+4(if“) 

3x  3y  3x  Re  3y 


0 

ay 


3U  +  9V  _  Q 


3x  3y 


P  Re(u  -51  +  v  -21)  =  EP(^)2  + 

ax  3y  3y  3  t 


where  the  dimensionless  parameters  are  defined 


Re*  =  Mk  (^0)2 


C  y 

p  =  P  s. 

K 


(2.8) 

(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 


U 

C„  T 
P 


(2.14) 


— 
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The  boundary  conditions 

—  —  _  \ 

u  =  1 ,  v  =  0  at  y  =  0 

u  =  0,  V  =  0  at  y  =  1  | 

\ 

p  =  0  at  x  =  0 

► 

p"  =  0  at  x  =  1 

/ 

T(0,y)  =  1  N 
T(x,0)  =  1  ) 

T(x,l)  =  1  ) 

2.3  TRANSFORMATION  OF  THE  MOMENTUM  EQUATION 

Let  the  stream  function  be  defined  by 


u 


dip 

dy 


v 


9X 


(2.15) 


(2.16) 


(2.17) 


(2.18) 

(2.19) 


where  after  SNYDER  [8] 
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*  =  K  l  f  (n)  6n  (2.20) 

n=o 

In  the  above  6  =  h/h ^  ,  n  =  y/h  ,  and  K  is  a  non  dimensi onali zing 
constant  to  be  specified  later. 

Expressing  the  appropriate  velocity  expressions  in  terms  of 
the  stream  function  gives 


K  v  -p'  * n-1 
u  ~  t—  >  f  6 

hi  n=0  n 


(2.21) 


3U 

3y 


T-  I  fn  6 
h^  n=0  n 


n-2 


(2.22) 


a  u  _  K  r  f 1 ' 1  .n-3 

3y^  =  hT  n=0  f|1  ^ 


(2„ 23) 


_au 

ax 


Kh 


^  •  _ 

r  _c  .n-2 


2  L  n 
h^  n=0  n 


f  6' 


Knh  y  f"  . 

~Z~  L  Tn  5 

h^  n=0  n 


n-2 


Kh 


l  n  fn  6 

h.  n=l  n 


n-2 


2 

'i 


(2.24) 


V  =  l  f'  60'1 

hi  n=0  n 


^  [  Bf  a-1-1 

hi  n-1  n 


(2.25) 


■ 


where  h  =  ~i—  -  constant, 
dx 

On  substituting  (2.21)  through  (2.25)  into  the  differential 
equation  (2.1)  the  following  is  obtained: 


2  ' 

ICh  p 

uK 


00 


I 

n=0 


f'  6n*2 
n 


.  K2h'p 
pK 


2  ' 
rh  p 

yK 


fn  4 


n-1 


+ 


9£ 

8X 


i  i  i 

f 

n 


Let  the  constant  K  be  such  that 


2  1 
rh  p 

yK 


which  yields 


K  = 


V 


h 


and  consequently  (2,26)  becomes 


■ 


16 


+ 


oo 

V 


L 


n=0 


f 

n 


n  f  6 
n 


n-2 

7 


oo 

y 

L 


n=0 


1 1 
f 

n 


6 


n-2  1 

I 


+ 


h3 

uK  3x 


6 


n-3 


(2.27) 


The  pressure  gradient  term  in  equation  (2.27)  is  also  assumed  to  have  the 
following  series  from  expansion  [8] 


h3  h3 

4H-4(c 


y  3X 


h' 


0 


C  -j  6 


C26‘ 


C363  + 


(2.28) 


where  C's  are  constants.  Combining  equation  (2.28)  and  (2.27)  a  set 
of  differential  equations  is  found  of  which  the  following  are  the  first 
four 


6 


6 


6 


6 


f0  +  Vo  K  C0 


+  f0Vvi  =!ci 


f2  +  2f0  f2  +  fl  fl  "  K  C2 


f3  "  f0f3  +  3f0  f3  flf2  +  2fl  f2  +  flf2 


(2.29) 

(2.30) 

(2.31) 


x 

K 


C 


3 


(2.32) 
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2.4  ASSOCIATED  BOUNDARY  CONDITIONS 

The  boundary  conditions  to  be  satisfied  by  the  momentum 
equation  are  (2.5)  to  (2.6). 

From  the  boundary  condition  (2,5),  a  set  of  conditions  for 
function  f  are  obtained  by  introducing  this  boundary  condition  into 
the  equations  (2 „ 2 1 )  and  (2.25).  These  are  the  conditions  that  are  to 
be  satisfied  by  the  differential  equations  (2.29)  to  (2.32). 


and 


where 


f0(0)  =  0  , 

f  (1)  =  — 

V  '  pK 

f,(0)  =  0  , 

f-,0)  =  0  , 


yo)  =  o 

I 

f0 (D  ■  0  , 

,  h.U  ^ 

V°>  -4- 

> 

yi)  =  o  , 


fn(0)  =  0  ,  fn(0)  =  0 
fpO)  -  0  ,  f>)  ■  o 


(2.33) 


(2.34) 


(2.35) 


n  2 j  3j  o.. 

Each  of  the  equations  (2.29)  through  (2.32)  for  the  f-functions 


is  of  third  order,  and  thus  three  boundary  conditions  will  be  required  for 


!  Q)  /  ,  0 


18 

obtaining  the  solution.  However,  each  equation  also  involves  an  un¬ 
known  constant  C  which  will  require  one  more  condition  on  function  f. 

This  necessary  condition  for  CQS  C-j  and  Cn  will  be  fg(l),  f-|(l)  and 
f  (1)  respectively.  The  C's,  which  depend  on  m  as  seen  from  equation 
(2.33)  will  determine  the  pressure  gradient.  The  pressure  then  will  be 
obtained  by  integration  of  expression  (2.28). 
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CHAPTER  III 

SOLUTION  TO  THE  GOVERNING  EQUATIONS 

3.1  NUMERICAL  SOLUTION  TO  THE  MOMENTUM  EQUATION 

The  transformation  of  the  momentum  equation  (2.1)  which  includes 
the  inertia  terms  for  the  slider  bearing,  yields  a  set  of  differential 
equations  (2.29)  to  (2.32)  along  with  the  boundary  conditions  (2.33)  to 
(2.35).  The  methods  of  solutions  to  these  differential  equations  with 
their  respective  boundary  conditions  are  described. 

(a)  Solution  to  the  first  term  in  stream  function  expansion  (fg) 

The  differential  equation  (2.29)  for  fg  is  of  third  order  and 
non  linear.  Its  non  linearity  made  it  difficult  to  obtain  a  closed 
form  solution.  It  was  therefore,  decided  to  apply  some  numerical 
technique  to  obtain  the  solution. 

The  differential  equation  (2.29)  and  the  boundary  conditions 
(2.33)  are  reducible  to  an  equivalent  initial  value  problem  involving 
the  solution  of  a  system  of  first  order  ordinary  differential  equations, 
subjected  to  prescribed  initial  conditions.  The  initial  conditions  to 
be  imposed  on  fg  and  fQ  are  obvious,  however,  the  initial  value  of 

I  I 

fg  is  not  readily  prescribed.  Rather  a  trial  and  error  method 
is  adopted. 

To  begin  the  computation  some  value  of  CQ  is  assumed.  Then  the 

1 1 

differential  equation  is  solved  with  the  various  assumed  values  of  f  (0) 
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I  I 

until  the  boundary  conditions  fg(0),  fQ  (0),  and  f  G(  3)  are  satisfied. 

I  I  I 

It  is  found  that  the  relation  between  fg  (0)  and  fQ  (1)  is  approximately 

I  I 

linear  which  results  in  a  very  good  approximation  of  f  (0)  after  only 

I  I 

two  iterations.  The  final  value  of  f  (0)  is  found  by  iteratively  cor¬ 
recting  its  previously  obtained  value. 

After  solving  the  differential  equation  for  the  assumed  Cq  the 
value  of  f q ( 1 )  is  compared  with  the  respective  boundary  conditions.  In 
case  of  disagreement  a  new  value  of  Cq  is  assumed  and  the  above  procedure 
is  repeated  until  the  boundary  conditions  are  satisfied.  It  was  found 
that  the  relation  between  Cq ' s  and  the  values  of  f q ( 1 )  obtained  in  the 
above  procedure  is  also  approximately  linear.  Therefore,  a  better  ap¬ 
proximation  of  Cq  was  made  possible.  The  final  value  of  Cq  is  found  by 
iteratively  correcting  its  previously  obtained  value. 

Gill's  variation  of  the  Runge-Kutta  fourth  order  method  is 
employed  as  an  integration  procedure.  The  subroutine  DRKGS  from  IBM 
System/360  scientific  subroutine  package  is  employed  to  obtain  the 
sol uti on . 

( b )  Solution  to  the  second  term  in  stream  function  expansion  (f-j J_ 

The  differential  equation  (2.30)  for  f^  is  of  third  order  and 
linear,  and  the  coefficients  involving  f q 1 s  has  already  been  obtained. 

The  solution  to  this  equation  with  the  boundary  conditions 
(2.34)  is  obtained  with  the  following  numerical  procedure. 

The  derivatives  in  the  differential  equation  are  approximated 


. 

no'rJjsupe  f,fitJri9'«$ftf,b  9rl}  nr  29vr^evn9b  9dl 
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by  the  finite  difference  approximations.  The  resulting  equation  can  be 
represented  by  the  recurrance  relation  as  suggested  by  Rodkiewicz  and 
Reshotko  [13] 


S  =  R  S  , ,  +  D 
n  n  n+1  n 


(3.1) 


where 


and 


-a 


n 


D  = 


R"  =  (bn  +  cnRn-l  +  dnVl  V2' 


e  -  (c  D  ,  +  d  R  0D  ,  +  d  D  0) 
n  v  n  n-1  n  n-2  n-1  n  n-2' 


n  <bn  +  cnRn-l  +  dnRn-l  Rn-2J 


(3.2) 


(3.3) 


where  a  ,  b  ,  c  ,  d  and  e  are  defined  in  Appendix  A. 
n  n  n  n  n 

For  n  =  1  we  have  from  the  boundary  condition  (2.34) 


R-|  -  S-|  -  D-|  -  0 


and 


h.U 


—  =  —  f-S  +4S  -S  )  =  — 
an  2H  1  H  ^2  V  K 

On  considering  equations  (3.4b)  and  (3.1)  for  n  =  2  we  obtain 


(3.4a) 


(3.4b) 


M 

R2  =  0.25  ,  D2  =  x  y 


(3.4c) 


For  n  >_  3  one  can  find  from  (3.2)  and  (3.3)  all  Rn  and  Dn 

respectively  as  we  will  have  Rn_2,  Rn-1’  Dn-2  anc*  Dn-T  Wlth  t*1ese 
values  one  can  obtain  from  (3.1)  variable  Sn  provided  we  have  S  +-|.  Let 
n  =  N  represent  the  location  for  which  the  upper  boundary  condition 

I 

prevails.  We  know  that  at  this  location  S  (1)  =  0  and  therefore 


f  ■  7H  <-3Sl+4S2-S3>  ■  0 


(3.5) 
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and  consequently 


SN-2  4SN-1  "  3SN 
Also  from  (3.1)  we  obtain 


(3.6) 


SN-1  ~  RN-1SN  +  °N-1 


(3.7) 


and 


SN-2  =  RN - 2SN - 1  +  °N-2 


Substituting  (3.6)  into  (3.8)  we  have 


3S  +  D 
^N  uN-2 


1-1  4  -  R 


N-2 


and  using  (3.9)  in  (3.7)  we  finally  obtain 


D, 


N-2 


-1  4  -  R, 


SN 


N-2 


4  -  R„  „  '  RN-1 


1-2 


(3.8) 


(3.9) 


(3.10) 


which  constitutes  the  starting  point  in  the  use  of  expression  (3.1). 

To  begin  the  computation,  some  value  of  the  unknown  constant 
C-|  in  the  differential  equation  (2.30)  is  assumed.  The  assumed  magnitude 
is  then  varied  until  all  the  boundary  conditions  are  satisfied.  These 
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computations  were  done  on  IBM  System/360. 

(c)  Solution  to  the  third  and  the  fourth  terms  in  stream  function 

expansion  (f^  and  fj 

The  differential  equations  (2,31)  and  (2-32)  are  of  third 
order  and  linear.  The  solutions  to  these  equations  with  their  respective 
boundary  conditions  are  obtained  by  using  the  numerical  procedure  as 
described  in  section  3. lb 0 

(d)  Pressure  distribution 

The  solution  of  differential  equations  (2.29)  to  (2.32)  yields 
the  values  of  Cq,  C-j  ,  ,  and  which  are  now  inserted  in  the  equation 

(2.28) o  Hence  giving  the  expression  for  the  pressure  gradient  along 
the  length  of  the  bearing.  This  expression  is  integrated  using  Simpson's 
Integration  Rule,  resulting  in  a  pressure  distribution.  Now  the  obtained 
value  of  pressure  at  the  outlet  is  compared  with  the  boundary  condition 
(2.6)  and  if  these  two  values  agree  within  the  specified  limit,  then  we 
have  obtained  the  solution  to  the  problem.  In  case  of  disagreement,  a 
new  assumption  of  the  mass  flow  is  made  and  the  above  procedure  is 
repeated. 

3.2  NUMERICAL  SOLUTION  TQ  THE  ENERGY  EQUATION 

The  solution  to  the  energy  equation  (2.11)  with  the  boundary 
conditions  (2.17)  is  obtained  by  the  method  of  finite  differences.  The 
derivatives  are  replaced  by  the  central  finite  difference  approximations, 
and  the  terms  are  arranged  to  obtain  (see  Appendix  B)  the  equation  in  the 
following  form. 


■ 

gnofs  Jrmbe'ie  9^u229"io  ed;  '10T  nc  229^q  j  9  pn.vrp  e-'neH  ,(8S,,S) 
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m  ,i 


(3.11) 


where  am  ,  b  ,  c  n  and  dm  are  defined  in  Appendix  B. 
m,n  m,n  m,n  m,n  rr 

The  boundary  conditions  become 


(3.12) 


(3.13) 


(3.14) 


In  order  to  obtain  the  solution  the  grid  presented  below  was 
employed.  The  spacing  interval  in  7  and  y  direction  was  taken  0.1  and 
0.05  respectively. 


FIG.  3.1  GRID  NETWORK 


■ 
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Equation  (3.11)  at  location  m  =  2  yields 


a2 ,n  T2 ,n-l  +  b2,n  T2,n  +  c2,n  T2,n+1  "  d2,n 

(3.15) 

and  the  boundary  conditions  become 

ii 

(3.16) 

ii 

(3.17) 

This  yields 

T2 ,2 9  T2,3* 
solution  of 

a  linear  system  of  N-2  equations  with  N-2  unknowns  namely 

.  T2  j^j  1  •  The  associated  coefficients  are  known  from  the 

the  momentum  equation  and  the  initial  condition  (3.14). 

a2 ,2T2 ,1  +  b2 ,2T2 ,2  +  C2,2T2,3  =  d2,2 

a2,3T2,2  +  b2,3T2,3  +  C2,3T2,4  =  d2,3 

(3.18) 

a2  ,N-2T2  ,N-3  +  b2  ,N-2T2  ,N-2  +  C2,N-2T2,N-1  =  d2,N-2 
a2  ,  N - 1 T  2  ,N-2  +  b2  ,N-1T2,N-1  +  C2,N-1T2,N  =  d2,N-l 


« 

11, ,  - 
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Restating  the  equation  (3/18)  in  the  matrix  form,  the  following 
is  obtained: 


b2 ,2  C2 ,2 

T2 ,2 

d2,2'a2,2 

a2,3  b2 ,3  C2  ,3 

T2  ,3 

d2 ,3 

e 

o 

T2,4 

d2  ,4 

a2,N-2  b2,N-2  C2  ,N-2 

T2,N-2 

d2,N-2 

a2  ,N- 1  b2  ,N - 1 

T2,N-1 

_ 

d2  ,N-l"C2 ,N-1 

The  above  matrix  is  a  tri diagonal  and  of  order  l\l-2.  Gauss's 
elimination  method  is  employed  to  obtain  the  solution. 

Once  T2  n  is  obtained  equation  (3.11)  can  be  solved  for 
m  =  3,4,5,...  using  the  same  procedure. 

The  subroutine  DGELB  from  IBM  System/360  scientific  subroutine 
packages  is  used  for  computations. 


, 


■ 
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CHAPTER  IV 

DISCUSSION  AND  CONCLUSIONS 


4.1  INERTIA  EFFECTS 

k 

The  momentum  equation  (2.1)  is  solved  for  Re  =  0,1,  0,25, 

0.5,  0.75  and  1.0.  For  comparison  the  computations  were  made  for  the 
case  of  neglected  inertia  terms.  The  dependent  variables,  pressure, 
load  capacity  and  shear  stress  at  the  moving  component  are  obtained. 

The  tabulated  results  and  plotted  graphs  are  given  in  Appendix  C. 

Tables  4.1  and  4.2  give  the  values  of  p"  and  W  at  different 
Re  with  and  without  inertia  terms  respectively.  Table  4.3  shows  the 
increase  in  load  capacity  due  to  the  presence  of  inertia  terms. 

Figures  4.1  and  4.2  shows  the  pressure  distribution  with  and 
without  the  inertia  terms  for  Re  =  0.25  and  1.0  respectively.  Figure 

k 

4.3  gives  the  pressure  computed  with  the  inertia  terms  for  Re  =  0.25, 
0.5,  0.75  and  1.0.  It  appears  from  these  curves  that  the  location  of 

k 

the  maximum  pressure  is  approximately  independent  of  Re  and  the  pressure 

* 

increases  as  the  Re  increases. 

Figure  4,4  shows  the  shear  stress  distribution  at  the  moving 

k 

component  for  Re  =  0.25,  0.5,  1.0.  The  shear  stress  increases  with 
the  increase  of  Re. 

In  Fig,  4.5  the  velocity  distribution  has  also  been  plotted 

k 

at  different  sections  along  the  bearing  for  Re  =  0,25,  0.5  and  1.0,  It 


' 


■ 
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is  noted  that  the  curves  of  the  dimensionless  velocity  become  more 

k 

concave  with  the  increase  of  Re. 

The  load  capacity  and  the  increase  in  load  capacity  due  to 

presence  of  inertia  terms  are  plotted  in  Fig.  4,6.  The  maximum  in- 

* 

crease  in  the  load  capacity  was  found  to  be  8,9%  at  Re  =  1,0. 

In  order  to  make  a  comparison  with  the  work  of  Kahlert  [6],  and 
Woodhead  and  Kettleborough  [9],  the  dimensionless  load  capacity  is 
redefined  according  to  [9],  The  pertinent  curves  are  plotted  in  Fig. 

4,7.  It  can  be  seen  that  the  present  results  are  in  very  close  agree¬ 
ment  with  those  of  Woodhead  and  Kettleborough. 

A  similar  computation  was  also  carried  out  by  considering 

only  the  first  two  terms  of  expression  (2.28).  Table  4,4  gives  the 

★ 

pressure  distribution  and  the  load  capacity  for  various  Re.  Fig.  4.8 
shows  the  curves  of  pressure  distributions  calculated  considering  the 
first  two  terms  only  as  well  as  considering  the  first  four  terms  of  the 
expression  (2.28)  for  Re  =  1.0.  Figure  4,9  shows  the  increase  in  load 
capacity  percentage-wise  due  to  the  additional  contribution  of  the  3rd 
and  the  4th  term  of  the  series.  It  is  noted  that  the  contribution  of 
the  last  two  terms  of  the  expression  (2,28)  is  approximately  74%  at 

k 

Re  =  1,0,  hence  their  effects  cannot  be  ignored  as  has  been  suggested 
by  SNYDER  [8].  Thus  it  would  seem  that  it  may  be  necessary  to  include 
at  least  four  terms  in  the  pressure  expression  (2,28)  to  obtain  more 


accurate  results. 
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4.2  CONVECTION  EFFECTS 

The  dimensionless  form  of  energy  equation  (2.11)  shows  that 
the  solution  of  the  differential  equation  depends  upon  three  important 
parameters  -  the  modified  Reynolds  number,  the  Prandtl  number  and  the 
Eckert  number.  For  computations  the  lubricant  with  P  =  122  is  selected 
here.  The  differential  equations  are  solved  with  and  without  convective 

•k 

terms  for  Re  =  0.25,  and  1.0.  The  corresponding  values  of  E  for 
Ts  =  1 50  F  are  0,000473  and  0.00757  respectively. 

Figures  4,10  and  4.11  show  the  temperature  distribution  in 
the  fluid  film  without  the  convective  terms  at  different  locations 

along  the  bearing.  The  maximum  temperature  rise  was  found  to  be  1 . 3°F 

- 

and  20,0~F  for  Re  =  0.25  and  1.0  respectively. 

The  temperature  distribution  in  the  fluid  film  with  the  con¬ 
vective  terms  is  plotted  in  Figs.  4,12  and  4,13.  The  maximum  temper- 

c  * 

ature  rise  was  found  to  be  0.4  and  2. OF  for  Re  =  0,25  and  1.0  respectively. 

Figures  4.14  and  4.15  show  the  temperature  gradient  with  and 
without  the  convective  terms  at  the  moving  and  stationary  components  for 
Re  =  0.25  and  1.0  respectively. 

The  results  i.e.,  the  temperature  distribution  and  the  temper¬ 
ature  gradient  obtained  with  and  without  the  convective  terms  for 
★ 

Re  =  1.0  indicates  that 

(i)  the  maximum  temperature  rise  in  the  film  is  approximately 
10  times  higher  for  the  case  of  neglected  convective  terms, 

(ii)  the  minimum  error,  due  to  neglecting  convective  terms,  in  the 


. 

. 

: '  : 

* 

, ag M  nr 


temperature  gradient  at  the  moving  and  stationary  surfaces  is  approxi¬ 
mately  1255%  and  1411%  respectively.  This  minimum  error  is  obtained 
by  the  following  expression 
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(— ) 

Without  Convective  Terms 


(— ) 

With  Convective  Terms 


(— ) 

'With  Convective  Terms 
ay 

Thus  it  would  seem  that  the  effect  of  convective  terms  in 

the  analysis  of  hydrodynamic  bearing  is  by  no  means  negligible  and 

★ 

their  omission,  particularly  at  high  Re  may  introduce  appreciable 
errors . 


4.3  TEMPERATURE  VARIATION  AT  HIGH  PRANDTL  NUMBERS 

The  momentum  equation  (2.8)  and  the  energy  equation  (2.11) 
are  also  solved  for  high  viscosity  lubricating  oils.  Figure  4.16 

■k 

shows  the  temperature  variation  for  P  =  2174.6.  The  Re  and  E  are 
0.04872  and  0.0082  respectively.  The  maximum  temperature  rise  in 
the  film  is  found  to  be  34°F. 

Furthermore  a  similar  computation  were  also  carried  out  for 
a  light  lubricating  oil,  P  =  347.8.  The  Re  and  E  being  0.305  and  .0082 
respectively.  The  temperature  variation  in  the  film  is  plotted  in 
Fig.  4.17.  The  maximum  temperature  rise  was  found  to  be  6.0°F. 

In  comparing  the  temperature  distribution  obtained  for  the 
above  two  cases  it  is  noted  that  the  lubricant  of  high  Prandtl  number 
develops  larger  temperature  variation  in  the  film. 
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4.4  CONCLUDING  REMARKS 

In  obtaining  the  solution  of  the  differential  equations  there 

were  certain  assumptions  made  vide  Chapter  II .  It  is  noted  from  the 

present  results  that  when  the  viscosity  of  the  lubricant  is  high,  and 

* 

particularly  if  Re  is  very  small,  large  temperature  variations  develop 
within  the  film.  For  such  cases  it  would  be  necessary  to  incorporate 
viscosity  dependence  on  temperature  in  the  differential  equations. 

This  would  bring  additional  difficulties  in  solving  the  differential 
equations  as  the  momentum  equations  will  be  now  coupled  with  the 
energy  equation  through  viscosity. 

A  further  assumption  made  was  that  the  boundaries  are  rigid. 
In  practise  this  may  not  be  true  for  the  case  of  bearings  developing 
very  high  pressures  and  temperatures .  For  such  cases  it  would  be 
necessary  to  include  the  equations  of  elasticity  for  the  boundaries 
in  addition  to  the  differential  equations  governing  the  fluid  flow. 
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APPENDIX  'A' 

RECURRENCE  RELATIONS  INVOLVING  THE  4-POINT 
COMPUTATIONAL  MOLECULE 

The  differential  equation  (2=30)  may  be  written  for  simplicity 
in  the  following  form 


S 


n 


+AnSn+BnSn  =G 


(A.  1 ) 


where  S  is  a  dummy  variable,  A  and  B  are  known  coefficients  and  G 

J  n  n 

is  a  constant. 

The  following  shape  of  the  computational  molecule  will  be 
used  in  the  numerical  method. 


n 

4 


H 


FIG.  A. 1 
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Let  H  be  the  step  size  in  n  di recti  on „  Then  [13]  we  may  write 


Sn  2H  ^Sn+1  “  Vl^ 


(A. 2) 


H' 


^Sn+1  "  3Sn  +  3Sn-l  '  Sn-2^ 


(A. 3) 


Substituting  (A. 2)  and  (A. 3)  into  the  equation  (A.1)  and  after  re¬ 
arranging  we  obtain 


d  +  T  Bn>  Sn+1  +  (HV3>Sn  +  <3  '  T  Bn>Vl 


'  Sn-2  ■  H  G 


(A. 4) 


The  expression  (A. 4)  can  also  be  written  in  the  following  general  form 


aS.n+bS  +  cS  n  +  dS  0  =  e 
n  n+1  n  n  n  n-1  n  n-2  n 


(A. 5) 


where 


an  ’  +  V  Bn> 


(A. 6) 


bn  ■  (H\  -  3) 


(A. 7) 


cn  ■  <3  -  T  Bn> 


(A. 8) 


d  =  - 
n 


(A. 9) 


1 . 0 
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(A. 10) 


Considering  the  boundary  conditions  in  conjunction  with  (A. 5)  a  recurrance 
formula  of  the  following  from  can  be  developed  [13] 


S 


n 


p  S 

'n  n+1 


+  D 

n 


. 
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APPENDIX  ' B 1 

ENERGY  EQUATION  IN  FINITE  DIFFERENCE  FORM 

The  energy  equation  (2.11)  is  approximated  using  a  4-point 
as  well  as  a  5-point  computational  molecule.  The  4-point  molecule 
approximation  is  only  employed  at  m  =  2  location,  and  for  m  =  3  onwards 
the  5-point  molecule  is  used  for  better  approximation. 

B.l  4-POINT  FINITE  DIFFERENCE  APPROXIMATION 

For  the  4-point  finite  difference  approximation  the  following 
computational  molecule  was  considered 


FIG.  B.l 
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Let  1  and  H  be  the  elementary  step  magnitude  in  x  and  y  direction 
respectively.  Then  [14]  we  may  write 


9_T 

dx 


T  (T  i  -  T  J 

I  m- 1  ,n  m,n 


(B.  1 ) 


il 

3y 


2H 


(B.2) 


1 

7 


+  T. 


m,n-l 


(B.3) 


Substituting  (B.l)  through  (B.3)  in  the  energy  equation  (2.11)  we  obtain 


U  * 

-  (1  +  2-  PRe 


vm,n^ 


m,n-l 


+  (2  + 


o  *  _ 

n  PRe  u 


1 


m,n_\  j 


m,n 


+  H  +  7  PRe*  Vn>  ,n+ 1  ■  h2[RE 

dy 


PRe  u  T 
+  _ m,n  m,n  j 


(B.4) 


which  could  be  written  as 


(B.5) 


where 


m,n 


-  O  +7pr1  Vn> 


(B.6) 
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(2  + 


* 


PRe  - 

-  Um  n 

m  ,n 


) 


(B.  7) 


n  vn) 


(B.8) 


3U  o 

(  m,ny4+ 

3y 


PRe  u 


m,n 


1 


T 

m,nj 


(B.9) 


B . 2  5-POINT  FINITE  DIFFERENCE  APPROXIMATION 

For  the  computational  molecules  indicated  below  in  Fig.  B.2 
we  may  wri  te 


FIG.  B.2 
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il 

dx 


21  ^3Tm+l,n 


+  T  ,  ) 

m-1  ,n' 


(B. 10) 


—  =  —  (T  -  T  ) 

-  2H  '  m+1 sn+l  'm+l,n-l' 


(B.ll) 


JlI  =  —  (t 

Sy2  H2  (  m+1>n+1 


2T  +  T  . ) 
m,n  m+ I ,n- I 


(Bo  12) 


Substituting  (B.10)  through  (B.12)  into  the  energy  equation  we  obtain 


PR*  Vl,n>  Tm+1  ,n-l  +  <2  + 


3H2PRe  u  .  -1  n 

m+  I  ,n \— 

21  '  m+1 ,n 


+  H  +JPR*  Vl.n>  Vl 


3U_ 


_  |,2rprfu  m+1  ,n \2 

>n+1  -  h  im-rz-) 


+ 


PRe  u 


m+1  ,n 


21 


(Bo  1 3) 


which  could  be  written  as 


am+l  ,nTm+l  ,n-l  +  bm+l  ,n^m+l  ,n  +  cm+l  ,nTm+l  ,n+l  ^m+1  ,n 


(B .  1 4) 


where 


m+1  ,n 


-  o  +  f  Vi,n) 


?  *  _ 
31-rPRe  u 


bm+l ,n  + 


m+1  ,n 


) 


(Bo  15) 


21 


(B .  16 ) 
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cm+l ,n 


(1  -  y  PRe  v  i  ) 
v  2  m+1  ,rr 


(Bo  1 7) 


9u 


m+1  ,n 


H2[pe  (  nTf^-n-)2  +  PRe 

3y 


m+l  ,n  (4T 


21 


m,n 


Tm-1  ,n^ 


(B . 18) 


APPENDIX  C 


TABLES  AND  CURVES  FOR 


CHAPTER  IV 


PRESSURE  DISTRIBUTION  WITH  INERTIA  TERMS 
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PRESSURE  DISTRIBUTION  WITHOUT  INERTIA  TERMS 
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TABLE  4.3 

INCREASE  IN  LOAD  CAPACITY 


* 

Re 

%  Increase  in  Load  Capacity 
Due  to  Presence  of  Inertia 
Terms 

0.1 

1.764 

0.25 

2.764 

0.5 

4.676 

0.75 

6.89 

1.0 

8.9 

Table  prepared  for  =  2 


PRESSURE  DISTRIBUTION  CONSIDERING  THE  FIRST 
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WITH  INERTIA  TERMS 


FIG.  4.1  PRESSURE  DISTRIBUTION 


WITH  INERTIA  TERMS 
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FIG.  4.2  PRESSURE  DISTRIBUTION 
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FIG.  4.3 


PRESSURE  DISTRIBUTION 
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FIG.  4.4  SHEAR  STRESS  DISTRIBUTION  AT  THE  MOVING  COMPONENT 
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FIG.  4.5  VELOCITY  DISTRIBUTION 
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FIG.  4.6  LOAD  CAPACITY 


INCREASE  IN  W  DUE  TO  THE  PRESENCE  OF  INERTIA  TERMS  -  PERCENT 
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FIG.  4.7  LOAD  CAPACITY 


INCREASE  IN  W  DUE  TO  PRESENCE  OF  INERTIA  TERMS  -  PERCENT 


SOLUTION  WITH  4 
TERMS  IN  (2.28) 
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FIG.  4.8  PRESSURE  DISTRIBUTION 
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FIG.  4.9  INCREASE  IN  LOAD  CAPACITY 


STATIONARY  COMPONENT 
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FIG.  4.10  TEMPERATURE  DISTRIBUTION  IN  FLUID  FILM 
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FIG.  4.11  TEMPERATURE  DISTRIBUTION  IN  FLUID  FILM 
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FIG.  4.14  TEMPERATURE  GRADIENT  ON  MOVING  AND  STATIONARY  COMPONENTS 
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FIG.  4.15  TEMPERATURE  GRADIENT  ON  MOVING  AND  STATIONARY  COMPONENT 


FIG.  4.16  TEMPERATURE  DISTRIBUTION  IN  FLUID  FILM 
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FIG.  4.17  TEMPERATURE  DISTRIBUTION  IN  FLUID  FILM 


